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SLOW  WAVE  CYCLOTRON  INSTABILITY 
IN  DIELECTRIC  LOADED  WAVEGUIDE  OF 
RECTANGULAR  CROSS  SECTION 


I.  INTRODUCTION 

The  instability  of  the  waveguide  modes  interacting  with  a  relativistic  electron  beam  guided 
by  an  applied  uniform  magnetic  held  B0  may  occur  due  to  two  different  mechanisms,  namely, 
cyclotron  maser  radiation1-2  (CMR)  and  Weibel  instability3.  Both  mechanisms  depend  on  an 
anisotropic  velocity  distribution  of  the  electrons.  The  instability  in  CMR  results  from  the 
azimuthal  bunching  of  electrons  due  to  the  relativistic  variation  of  their  masses.  Weibel  insta¬ 
bility,  on  the  other  hand,  is  nonrelativistic  in  nature  and  originates  from  axial  bunching  of  the 
electrons  induced  by  the  Lorentz  force  VL  x  B(U  where  is  the  velocity  of  electrons  perpen¬ 
dicular  to  B0  and  fl(U  the  wave  magnetic  field.  The  phase  bunching  tends  to  be  destroyed  if 
there  is  any  spread  in  the  velocity  of  the  electrons.  In  a  recent  study4,  it  was  shown  that  the 
two  mechanisms  always  combine  to  offset  one  another.  The  azimuthal  bunching  predominates 
for  waves  with  phase  velocity  (v^)  greater  than  the  speed  of  light  in  vacuum  (c)  while  the 
axial  bunching  predominates  if  v^  <  c.  Thus  cyclotron  maser  radiation  (CMR)  is  generated  in 
fast  wave  modes  and  the  Weibel  instability  in  slow  wave  modes.  Weibel  instability  is  capable  of 
producing  amplification  over  a  wide  bandwidth.  In  this  paper  we  consider  slow  wave  cyclotron 
instability  in  a  dielectric  loaded  waveguide  of  rectangular  cross-section.  In  this  geometry,  the 
instability  can  be  studied  in  the  fundamental  mode  which  couples  efficiently  with  external  cir¬ 
cuit. 
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In  Section  II  we  derive  the  modified  dispersion  relation  of  the  dielectric-loaded  waveguide 
modes  in  the  presence  of  a  weak  electron  beam.  The  response  of  the  electrons  is  determined 
by  solving  linearized  Vlasov  Equation.  Effect  of  a  spread  in  VL  and  Vt  (components  of  electron 
velocity  perpendicular  and  parallel  to  the  applied  magnetic  field  Bo)  is  also  considered.  We 
assume  that  the  spread  in  VL  and  Vt  arise  from  a  spread  in  the  angle  at  which  the  electrons 
enter  the  magnetic  field  Bo  but  all  the  electrons  have  the  same  energy.  Results  are  discussed  in 
Sec.  III. 

II.  MODEL  AND  FORMALISM 

Figure  1  shows  an  end-view  of  the  configuration  under  consideration.  An  annular  beam 
of  electrons  uniformly  distributed  over  a  cylinder  of  radius  r0  is  guided  by  a  uniform  magnetic 
field  B0  along  helical  trajectories  inside  a  waveguide  partially  filled  with  a  dielectric.  The  origin 
of  the  coordinate  system  is  at  the  center  of  the  guide.  Wave  propagation  takes  place  along  z- 
axis.  The  field  B0  is  parallel  to  the  z-axis.  Region  I  (-  d  <  x  <  d)  is  vacuum  and  region  II  (- 
Ly  <  x  <  —  d,  d  <  x  <  L\)  dielectric.  The  linear  response  of  the  electrons  is  obtained  under 
the  following  simplifying  assumption:  (i)  the  beam  is  sufficiently  tenuous  so  that  its  self  fields 
may  be  neglected,  (ii)  the  rf  fields  are  first  order  with  respect  to  B0  and  the  perturbed  electron 
distribution  function  /*"  (7,p,t)  caused  by  these  fields  is  of  first  order  with  respect  to  the  equili¬ 
brium  distribution  function  f0  (7,p),  and  (iii)  the  distribution  function  and  the  rf  fields  are 
independent  of  y-coordinated. 

The  modes  of  the  dielectric  loaded  waveguide  cannot  be  classified  as  TE  or  TM  with 
respect  to  the  propagation  direction.  However,  they  may  be  classified  as  TEx  or  TMx  with 
respect  to  the  x-axiss.  These  modes  are  also  classified  as  logitudinal  section  electric  (LSE)  and 
longitudinal  section  magnetic  (LSM)  modes,  respectively.  There  is  no  TMx  mode  which  is 
independent  of  the  y-coordinate.  The  modes  may  be  divided  into  odd  and  even  symmetries 
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with  respect  to  the  reflection  at  yz-plane.  It  will  show  that  the  even  and  odd  modes  couple, 
respectively,  to  the  odd  and  even  order  cyclotron  modes.  Since  our  primary  concern  is  coupling 
to  the  fundamental  cyclotron  mode,  we  give  detailed  calculation  for  the  waveguide  modes  of 
even  symmetry  and  mention  the  results  in  the  case  of  modes  of  odd  symmetry.  For  L2  <  2Lt, 
the  fundamental  waveguide  mode  belongs  to  even  symmetry.  In  accordance  with  assumption 
(iii)  we  restrict  our  attention  to  TExM  modes.  These  modes  with  even  symmetry  are  given  by* 
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and  the  unperturbed  dispersion  relation  is 

t*n  li*d  —  In,,  col  lnnt.  (6) 

n  and  €  are  respectively  the  permeability  and  dielectric  constant  of  medium  II.  The 

corresponding  vacuum  values  are  po  and  €0-  Slow  waves  exist  for  mo<o«  <  k,  <  pew.  In  this 
case  is  imaginary.  We  wilt  assume  that  m  “  M o-  A  typical  dispersion  curve  and  field  distri¬ 
bution  for  the  fundamental  mode  are  shown  in  Fig.  2a  and  2b,  respectively.  A  bend  in  the 
dispersion  curve  occurs  at  k,  -  k,  such  that  w  “  ck,.  From  Eq.  (4)  and  (6),  for  the  funda¬ 
mental  mode 


w/2r(€/€„-  !)•'*. 


(7) 
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For  kj  <  k,  we  have  u  >  ck  (fast  wave)  and  for  k,  >  k„  u  <  ck  (slow  wave).  For  fast  wave, 
the  field  is  maximum  at  the  center  of  the  wave  guide  while  for  slow  waves  the  maximum 
occurs  in  the  dielectric. 


In  the  presence  of  a  weak  electron  beam,  the  results  given  in  Eqs.  (l)-(7)  will  be  slightly 
modified.  Maxwell’s  equations  yield  for  B„ 


+  €(*)  k> 


!<•> 
r.  On 


ifikjyiu 
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where  j,u  is  the  perturbed  current  density  and 


«(x)  -  *o.  0  <  |x|  <  d 

-  «,  d  <  |x|  <  Lx.  (9) 

B$m  may  be  expanded  in  terma  of  the  unperturbed  functions  0Hjru0,m.  Multiplying  Eq.  (8)  by 

oB,%,  integrating  over  x  and  y  and  neglecting  the  off-diagonal  terms  on  the  left  hand  side,  we 

obtain 
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where  the  dimensionless  quantity  0,  is 
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In  writing  the  limits  of  the  x  -  intergration  on  the  right  hand  side  of  Eq.  (10),  we  use  the  rela¬ 
tion  that  Jy(u  —  0  for  |x|  >  d.  The  perturbed  current  j,u  is  given  by 

jym  ~-e  f  Vyfu(r,fit)dp 

~-‘fo  5* und>f"(7.p,i)dp.  (12) 


The  perturbed  distribution  function  /<n(7,J»lf)  is  obtained  from  linearized  Vlasor  equation  by 
integrating  along  the  unperturbed  orbits.  The  result  is  standard7  and  can  be  written  as 


/"(**/)  - 1  J'm  dt' t U)ty.n  +  .  vr/0(7\r). 


(13) 
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where  the  integration  is  along  the  unperturbed  orbits.  The  primed  quantities  7'  and  p '  are 
treated  as  functions  of  t'  while  7 and  p  are,  respectively,  the  values  of  7'  and  p'  at  /'  ™  t.  In 
Eq.  (13)  we  need  to  use  only  the  lowest  order  fields  given  by  Eqs.  (l)-(3).  We  now  assume  a 
specific  form  for  the  equilibrium  distribution  function  fo(7,p)  constructed  from  the  constants  of 

motion,  Pl,p,,x  —  x  +  Pl—  sin  <t>  and  y  —  y  +  — cos  (x  and  y  locate  the  guiding 
m0iio  niqIi  o 


go(PL-P^'  (l4) 


center  of  an  electron  and  ft0  “  eB/m<): 

N  (g(r0  -  x)  -  9(-r0  -  x)}  (g(r0  -  y)  -  »(-r0  -  j>)| 

2tr  r0  Vro  “  x2 

where  J  g0  (  pL,  pL)  dp  -  1  This  represents  an  annular  beam  of  electrons  uniformly  distri¬ 
buted  on  a  cylinder  of  radius  r0.  N  is  the  number  of  electrons  per  unit  length  of  the 
waveguide.  Since  we  consider  waveguide  modes  which  are  independent  of  y,  there  is  no  loss  in 
generalization  in  assuming  that  y  is  distributed  uniformly  in  the  range  -  r0  to  r0. 


The  modified  dispersion  relation  for  the  waveguide  in  the  presence  of  the  electron  beam 
may  now  be  obtained  from  Eqs.  (10)-(14).  The  methods  for  evaluating  the  integrals  in  these 
equations  are  standard7  and  the  result  may  be  written  in  the  form 


wl  2  _  2|  4rrMl  +  J0(2/,nr0)i 
c2  **J"  2LiL2B„ 
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l«-00  y 


(w  -  *2v2) 
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- X 
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S«o 
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1 

where 

y  -  (1  +  pf/mfic2  +  p22/m<?c2)l/: 
v  -  Ne2/Air  €0  m0c2. 
flo  “  eBg/ Wq, 

x  -  li„pL/m0(lo.  d6) 

and  J,(x)  the  Bessel  function  of  order  s. 
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We  now  specialized  to  a  particular  distribution  in  velocity  space  given  by 

go(PL-Pz)  “  /2<4'')  8(po  -  pz)6(p$  -  Pi  ~  ft1).  07) 

where  pi  “  p20  +  fto  and  the  normalizing  constant  AT  is  such  that  f  g0  dp  -  1 .  Thus 

Equation  (17)  represents  a  monoenergetic  beam  with  a  spread  in  the  angle  at  which  the  elec¬ 
trons  enter  the  field  B0 .  From  Eqs.  (IS)  and  (17)  we  then  obtain 


(“  "  l'n  "  " 


4 vK  f’o.  «  (1  ~  (-!)») 

LtLiyoT'Jo  P:  2 
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Oo  . 


where 


"  w(l  +  Jo  (2l/„ro)} 

X  -  Z/aVft?  ~  Pz2/m0n  o-  (20) 

For  waveguide  modes  of  odd  symmetry  an  expression  similar  to  Eq.  (19)  is  obtained  except 

that  the  factor  (1  -  (-  l)5)  becomes  {1  +  (-  l)5)  and  in  the  expression  for  S„  in  Eq.  (11)  sin 

2lt„d  changes  to  -sin2 lM  and  sin 2l,„d  to  cos2l,„d.  Thus  even  symmetry  modes  couples  with 

cyclotron  modes  having  odd  values  of  s  and  vice  versa.  The  unperturbed  dispersion  relation  in 

Eq.  (6)  is  replaced  by  4,  cos  lt„d  -  - 1,,„  cos  lUnt  for  odd  symmetry  modes. 


The  length  L]  can  be  scaled  out  of  the  dispersion  relation  in  Eq.  (19)  by  introducing  the 


following  normalization  scheme: 
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(0  length  normalized  to  L\  (e.g  r0  -  r0/Lt) 

(ii)  frequency  normalized  to  c/Lt  (e.g.  5  -  wLjc) 

(iii)  momentum  normalized  to  m0c  (e.g.  Pi  -  pjm0c). 

Other  quantities  such  as  kz.  and  V  are  to  be  normalized  consistently  with  the  proceeding 
procedures  (e.g.,  kz  »  kzL 7„„  —  lthL u  v,  —  v,/c,  vx  -  \Jc).  In  terms  of  these  dimension¬ 
less  variables,  the  dispersion  relation  (Eq.  (19))  may  be  written  as 


where  x  —  l/KyJpo  -  pzl ft 0-  For  slow  wave  7/„  (and  x)  is  imaginary.  In  this  case  J,(x)  - 
y,(/|x|)  —  //j(|x|)  where  Is  is  the  modified  Bessel  function.  Also  sin  li„d  in  ihe  expression  for 
0.  changes  to  isinh\Ji„\d. 


where  x0  -  7Mpi0/n0. 
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For  slow  waves  <■>  <  kz  and  the  second  term  inside  the  square  bracket  in  Eq.  (22)  is  posi¬ 
tive.  Instability  can  occur  if  the  first  turn  is  negative  i.e.  w  <  s^-  +  k2  v2  since  v,  <  <H/kz. 

y  o 

For  cyclotron  maser  radiation,  however,  u  >  k2  and  the  second  turn  drives  the  instability  and 

the  condition  5  >  s — —  +  kz\2  is  required  for  azimuthal  bunching. 

y  o 

III.  RESULTS  AND  DISCUSSION 

The  dispersion  relation  Eq.  (21)  is  to  he  analyzed  numerically.  For  a  travelling  wave 
amplifier  it  is  convenient  to  treat  <5  as  real  and  kz  complex.  The  spatial  growth  of  the  rf  field  is 
then  given  by  the  gain  coefficient  ( g ): 

gUn  dB / cm)  -  -  8.69  Im  (fcr)/L,  (23) 

where  L,  is  in  cm.  Gain  occurs  if  Im  kt  <  0. 

In  Equation  (21)  there  are  ten  free  parameters  v,  L2,  rf,  70,  e/«0,  y0(  vi0  /v.0,  ft0, 
»  and  s.  We  choose  the  lowest  modes  with  n  -  1  and  s  —  1.  It  is  found  that  the 
bandwidth  of  the  amplifier  increases  with  decrease  in  d  and  increase  in  «,  =  t/t0.  If  rf  is  made 
too  small,  however,  it  is  difficult  to  accommodate  the  electron  orbits  is  region  /  of  Fig.  1.  For 
efficient  interaction  between  the  electrons  and  the  rf  field,  v2  should  be  nearly  equal  to  the 
group  velocity  of  the  waves.  vf  decreases  with  increase  in  «.  Too  low  a  value  of  v2  (and  vx  is 
not  desirable.  For  numerical  calculation  we  choose 

rf  -  .6,  «r  -  20,  70  -  .1  and  L2  -  2d. 

Figure  3  shows  a  plot  of  Im  kt  as  a  function  of  5  for  different  values  of  \p./p:(l  at  y0  ~  10617 
(v20“  -2,  vi0/v20  “  1.35)  and  ftjyo  "  -72.  The  bandwidth  and  gain  decrease  as  the  velocity 
spread  increases.  Each  value  of  Ap2  /pj0,  the  magnetic  field  (i.e.,  fl0)  has  to  be  optimized  for 
maximum  bandwidth.  In  Fig.  4  we  show  Im  0cz)  as  a  function  of  Z>  for  different  values  of 
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&P:/P:o  with  the  corresponding  fl0  for  maximum  bandwidth.  The  values  of  Ap./p^,  fl0  and 
the  bandwidth  AJo!  <3  are  given  in  Table  I.  Fig.  5  shows  a  plot  of  AJo  /Jo  vs  Ap,/p:0 ■  The 
bandwidth  decreases  rapidly  from  64%  at  Apz/pz0  =  0  to  12%  at  Ap./pz0  =  .05.  The  calculated 
bandwidths  correspond  to  a  change  in  gain  coefficient  given  by  the  relation  [gmax  -  g(c5)]/gma), 

=  In  Figs.  6-8,  we  plot  Im  k ,  as  a  function  of  Jo  to  show  the  changes  in  gain  and 

21og(,10 

bandwidth  with  variation  in  the  parameters  ft0,  v.0  and  the  ratio  a  —  vi0/v.0,  respectively.  It 
is  seen  that  the  gain  and  the  bandwidth  are  sensitive  functions  of  fl0  and  vj0  and  increases 
slowly  with  increase  in  a.  The  ratio  a,  however,  cannot  be  made  too  large  since  for  a  given  v.0 
and  ft0,  the  Larmor  radius  rL  increases  with  a. 

Table  I 


bPs/PiO 

flo/vo 

A  to/to 

fL 

0 

.680 

.64 

.3970 

.001 

.705 

.45 

.3830 

.01 

.720 

.28 

.3750 

.02 

.728 

.205 

.3711 

.03 

.734 

.165 

.3678 

.04 

.737 

.140 

.3664 

.05 

.739 

.125 

.3654 

Zj/2  -  7i 

v.0  -  2. 

’  6-  «, 
v10/v.0- 

-  20, 

1.35,  r 

Jo  -  •  1  • 

-  .002. 

>0  -  1  062 

To  convert  normalized  quantities  fin  and  to  to  design  parameters,  the  following  relations 


should  be  used.  For  a  fixed  frequency  /  (in  Hz) 


,  <oc  „ 

L^u?cm’ 


Ha 


n„ 

1.707  ——  kOe, 
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while  for  fixed  magnetic  field  H0  (in  kOe) 

L{—  1.707  x  Cl<JH0cm 
f  —  md2itL\  Hz. 

We  give  below  the  design  parameters  of  a  slow  wave  amplifier  for  operation  at  frequency  /  — 
24  GHz.  For  v  -  .002,  \l0  -  .2  and  vi0/v,0  =  1.35,  the  parameters  are:  beam  energy  of 
F„(kV)  —  31.52,  beam  current  Ib  (amp)  =  3.41,  2Z.(  -  .426  cm,  2d  —  L2  *  -256  cm,  r0  — 
.0213  cm,  H0  —  6.245  k  0  e,  e/e0  —  20.  In  this  case  Aw/w  —  .165  and  gmat  “  1.346  dB/cni 
when  Apj/p.o  —  .03. 

IV.  CONCLUSION 

We  have  shown  that  slow  wave  cyclotron  instability  may  be  used  to  obtain  broadband 
amplification  in  the  millimeter  range.  The  bandwidth  deteriorates  with  increase  in  the  velocity 
spread  of  the  electrons.  The  velocity  spread  should  not  exceed  2%  to  obtain  a  bandwidth  of 

30%  or  more.  Our  results  are  based  on  a  linear  theory.  A  nonlinear  analysis  is  necessary  to 
calculate  the  efficiency  and  the  saturated  output  power  level  of  the  system. 
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F*.  6  -  Gain  vs  frequency  for  three  different  values  of  II,  at  Afc/?»  ”  0.2.  Other  parameters  are  same  as  in 
Fig.  3.  Curve  I:  njy„  -  .71  and  A5/5  -  .16;  curve  2:  Hjy,  -  .726  and  Aw/5  -  .2;  and  curve  3:  Hjy0  - 
.74  and  Aw/fi  “  .15. 
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a 

a 

a 

a 


•  1.45,  AS/S*.  216 

■  1.40,  A  w/w-  .211 

■1.35,  .205 

■  1.30,  A«5/£«  .200 


1  »  1  »  I  «  1  j _ I  »  I 

.0  l.l  _  1.2  1.3  1.4  1.5 

w 

lm  i,  v»  5  for  different  vehiei  of  o  el  lijy,  ~  728 
Other  p«r»  me  ten  ere  Mine  et  in  Fit.  3. 


